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Abstract 

The  impetus  of  this  paper  is  to  analyse  numerically  and  experimentally  the  flow  and  heat 
transfer  characteristics  of  buoyancy-driven  air  convection  behind  photovoltaic  panels.  Both 
convection  and  radiation  heat  exchanges  are  considered  as  the  heat  transfer  mechanisms  by 
which  the  thermal  energy  is  transferred  into  the  air.  Numerical  and  experimental  results  are 
obtained  for  a  channel  of  7.0  m  in  height  and  the  channel  walls  are  separated  by  a  distance  of 
0.23  m.  In  the  experiment  heat  is  supplied  to  the  air  gap  from  heating  foil  attached  to  one  of 
the  vertical  walls.  Different  input  heat  fluxes  and  emissivity  of  the  bounding  surfaces  are 
considered  in  order  to  show  their  effect  on  the  heat  transfer  across  the  air  layer.  Detailed  studies 
of  the  flow  and  thermal  fields  in  the  air  are  presented  in  order  to  explore  the  thermal  behaviour 
of  the  air  in  the  channel.  Velocity  and  temperature  profiles  of  the  outlet  air  and  the  surface 
temperature  of  the  heated  and  insulated  wall  is  presented.  The  numerical  results  agreed  well 
with  the  experimental  measurements.  ©  1998  Elsevier  Science  Ltd.  All  rights  reserved. 

Key  words:  Photovoltaic  panels;  Turbulent  flow;  Heat  transfer;  Vertical  channel;  Natural  convection; 
Radiation 


Nomenclature 

A  area 

B  specific  buoyancy  flux 
Cp  specific  heat 
d  wall-to-wall  distance 
F  shape  factor 
g  acceleration  due  to  gravity 
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Gk  buoyant  production  of  k 
h  convective  heat  transfer 
H  height 

k  turbulent  kinetic  energy 
k  friction  factor,  see  eqn  (13) 
n  surface  unit  normal  vector 
q  total  amount  of  heat 
q"  heat  flux 

q"n  uniform  input  heat  flux 
p  pressure 

Pk  shear  production  of  k 
Pr  Prandtl  number 
Ra  Rayleigh  number 
T  temperature 
w  width  of  the  air  gap 

u,  v  velocity  components  in  the  v-,  y-direction 

U  flow  average  velocity 

x,  y  co-ordinates  in  the  Cartesian  system 

Greek  symbols 
a  configuration  factor 
aH  thermal  diffusivity 

thermal  expansion  coefficient 
d  Kronecker  delta 
e  emissivity 

£  turbulent  energy  dissipation  rate,  see  eqn  (8) 
a  thermal  conductivity 
X  friction  factor 
p  dynamic  viscosity 
v  kinematic  viscosity 
p  density 

o  Stefan-Boltzmann  constant 

Subscripts 
a  ambient 
c  insulated  wall 
conv  convection 
entr  entrance 
h  heated  wall 
i  inlet 
o  outlet 
rad  radiation 
t  turbulent 

1.  Introduction 

The  solar  cell  operation  is  based  on  the  ability  of  semiconductors  to  convert  sunlight 
into  electricity  by  exploiting  the  photovoltaic  effect.  There  is  a  growing  interest  in 
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studying  heat  and  mass  flow  in  air  gaps  behind  Photovoltaic  panels  (PV).  One  impetus 
for  this  type  of  work  is  the  interest  in  hybrid  systems  (i.e.  the  combined  generation  of 
heat  and  electricity).  A  hybrid  system  consists  of  a  PV  panel  with  a  cavity  (air  gap) 
between  the  PV  panel  and  the  building  fabric.  This  air  gap  is  a  natural  draft  ‘chimney’ 
and  the  flow  in  the  air  gap  is  driven  by  buoyancy  (heat)  and  the  wind-induced  pressure 
difference  between  the  top  and  bottom  of  this  chimney.  The  PV  panel  is  heated  by 
the  incident  solar  radiation  and  from  the  PV  panel  heat  is  transferred  to  the  air  gap 
by  convection  and  radiation.  Radiative  heat  transfer  carries  energy  across  the  air  gap. 
The  net  radiation,  absorbed  by  the  otherwise  unheated  wall,  is  in  turn,  transferred  to 
the  fluid  by  convection.  Thus  radiation  activates  convection  on  the  otherwise  unheated 
surface. 

Another  impetus  is  the  fact  that  by  reducing  the  solar  cell  temperature  the  efficiency 
of  the  cells  will  increase.  The  efficiency  of  some  types  of  solar  cells  is  very  much 
dependent  on  the  operating  temperature  of  the  cell.  For  example  the  relative  tem¬ 
perature  coefficient  of  crystalline  silicon  solar  modules  is  in  the  range  0.4-0. 6%  K. 
With  13%  absolute  conversion  efficiency  this  corresponds  to  an  absolute  temperature 
coefficient  between  0.031  and  0.046%  K.  Therefore  a  reduction  by  20°C  will  give  an 
increase  in  efficiency  between  0.6  and  1%. 

As  a  result,  knowledge  enlightening  the  actual  physical  conditions  in  air  gaps  in 
ventilated  facades  is  crucial  for  design  engineers.  The  measurement  conditions  are 
very  difficult  due  to  the  presence  of  strong  temperature  gradients  and  the  exchange  of 
heat  between  surfaces  by  radiation.  At  the  same  time  the  strong  heat  transfer  mech¬ 
anisms  pose  problems  when  using  Computational  Fluid  Dynamics,  CFD.  In  real 
buildings  there  are  almost  transient  conditions  and  both  the  heat  input  to  the  air  and 
the  flow  rate  in  the  air  gap  vary  with  time.  The  instantaneous  heat  input  to  the  air 
can  hardly  be  measured.  The  approach  taken  in  this  study  is  to  combine  different 
techniques,  CFD  and  experiments  which  shall  be  seen  as  complementary  techniques. 
The  basic  approach  is  to  seek  solutions  for  simplified  boundary  conditions.  For 
example,  only  stationary  conditions  are  explored  and  the  heat  input  to  the  air  is 
controlled  and  known.  In  the  CFD  predictions  the  heat  input  is  given  as  a  boundary 
condition,  while  in  the  experiments  the  heat  input  is  known  by  using  a  heating  foil 
and  a  well  insulated  structure.  The  solutions  for  simple  boundary  conditions  are 
asymptotic  solutions  to  the  more  general  transient  cases  occurring  in  real  buildings 
and  they  will  shed  light  on  the  range  of  conditions  occurring  in  ventilated  facades. 

The  techniques  are  complementary  in  the  sense  that  for  example  detailed  tem¬ 
perature  profiles  can  be  predicted  by  CFD  but  are  difficult  to  measure.  On  the  other 
hand,  experiments  can  cope  with  transitional  flows  which  CFD  cannot. 


2.  Description  of  physical  model 

In  studies  of  air  flow  and  heat  transfer  between  solar  cells  and  the  air  gap  in  hybrid 
systems  one  can  to  some  extent  mimic  the  properties  of  solar  panels  by  inserting  a 
heating  foil  in  an  insulated  channel.  By  choosing  a  thick  insulation  one  can  ensure 
that  most  of  the  heat  generated  by  the  heating  foil  is  transferred  to  the  air  in  the 
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channel.  This  approach  enables  a  systematic  study  of  the  relation  between  input  heat 
and  various  flow  and  heat  transfer  properties.  A  vertical  channel  with  an  air  gap  with 
dimensions,  w  =  1.54  m,  d  =  0.23  m,  and  height,  H  =  6.5  m  is  built  in  the  Laboratory 
for  Ventilation  and  Air  Quality,  at  the  Royal  Institute  of  Technology,  Gavle,  Sweden. 
The  channel  is  insulated  with  0.29  m  mineral  wool.  To  minimise  end  effects  the  width 
of  the  heating  foil  is  larger  than  the  width  of  the  air  gap.  The  heating  foil  is  460  mm 
wide  and  0.5  mm  thick.  The  heating  elements  in  the  foil  consist  of  electrically  con¬ 
ducting  transverse  stripes  8  mm  wide  spaced  2  mm  apart.  The  heating  foils  are 
attached  with  double  adhesive  tape.  The  heating  foil  covers  the  whole  space  of  the 
walls.  The  heating  foil  is  powered  by  a  variable  transformer.  The  electrical  power 
supplied  to  each  section  is  measured  by  means  of  a  power  transmitter  which  converts 
the  electrical  power  to  a  voltage  simply  presented  on  a  voltmeter.  The  relative  error 
in  recorded  power,  Pm,  is  assessed  to  be  less  than  2%.  The  tests  are  carried  out  for 
the  heat  fluxes,  20,  50,  100,  200  and  300  W/m2,  respectively. 

The  surface  temperature  on  the  walls  is  measured  by  means  of  foil  thermocouples 
attached  with  double  adhesive  tape.  Measuring  of  the  air  temperature  between  the 
walls  is  done  by  using  small  thermocouples  made  of  0.2  mm  diameter  wires.  The  error 
is  0.3°C  absolute,  and  the  relative  error  between  measuring  points  is  estimates  to  be 
less  than  0.1  °C. 

The  flow  rate  in  the  channel  is  measured  by  using  tracer  gas  with  the  so  called 
constant  flow  technique.  Gas  is  introduced  at  the  bottom  of  the  channel  with  a 
known  volume  flow  rate  and  the  volume  concentration  at  the  exit  is  recorded.  The 
concentration  is  read  with  an  infrared  gas  analyser  giving  the  volume  concentration 
in  parts  per  million  [ppm].  Tracer  gas  technique  is  robust  and  accurate  technique.  A 
review  of  tracer  gas  techniques  with  a  statement  of  the  accuracy  is  given  in  Etheridge 
and  Sandberg  [1]. 


3.  Description  of  numerical  model 

The  geometry  of  the  air  gap  and  the  relevant  parameters,  as  well  as  the  boundary 
conditions  considered  here  are  shown  in  Fig.  1.  The  origin  of  the  rectangular  Cartesian 
co-ordinate  system  (x,y)  is  located  on  the  bottom  line  of  the  channel  entrance. 

Vertical  surfaces  represent  a  heated  wall  with  a  height  of  6.5  m  exposed  to  a 
constant  uniform  input  heat,  q"n,  and  an  insulated  7.0  m  high  wall.  The  width  of  the 
walls  is  assumed  to  be  large  in  comparison  to  other  dimensions,  thus  the  problem  can 
be  considered  two-dimensional.  The  bottom  surface  is  also  well  insulated  and  assumed 
to  be  a  blackbody.  The  air  enters  at  a  temperature  equal  to  the  ambient  temperature, 
T.d,  at  the  lower  left  hand  side  and  the  buoyancy-induced  flow  passes  over  the  heated 
and  insulated  walls  and  exits  at  the  top  of  the  gap.  These  are  the  realistic  thermal 
boundary  conditions  which  can  be  considered  for  the  case  presented  in  this  study. 

The  flow  is  assumed  to  be  steady  with  negligible  viscous  dissipation.  The  thermo¬ 
physical  properties  of  the  air  are  assumed  to  be  constant,  except  for  the  buoyancy 
term  of  the  momentum  equations,  i.e.  the  Boussinesq  approximation,  see  Gray  et 
al.[2]  The  radiation  surfaces  are  assumed  to  be  grey-diffuse  and  the  air  does  not 


B.  Moshfegh,  M.  Sandberg / Renewable  and  Sustainable  Energy  Reviews  2  ( 1998 )  287-301 


291 


Insulated  wall 
dT/dx-0 


Fig.  1 .  Geometry  and  boundary  conditions  for  the  numerical  model. 


interact  directly  with  the  radiative  process.  On  such  modelling  assumptions,  the 
governing  equations  for  conservation  of  mass,  momentum,  energy  and  turbulent 
kinetic  energy  as  well  as  turbulent  energy  dissipation  are  obtained  as  follows: 


Equation  of  the  conservation  of  mass 

V-ii  =  0  (1) 

Equation  of  the  conservation  of  momentum 

pV  •  uu  =  —  dp/dx  +  V  •  (p  +  pt)Vw  +  V  •  (p  +  iit)(du/dx)  (2) 

pV-uu  =  -  dpldy—gpfi{T—  T.d)  +  V  •  (p  +  pt)Vr  +  V  *  (p  +  gt)(du/dy)  (3) 

Equation  of  the  conservation  of  energy 

pCpV-wE=  V-(2  +  2t)VE  (4) 


where  the  heat  transfer  condition  on  each  segment  of  the  boundary  of  the  fluid  and 
the  heat  exchanged  between  radiating  boundaries  regardless  of  the  fluid  is  calculated 
by: 
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q"  =  -(AVT)n(s)  =  q'&s)  +  £/",nvCv)  +  qmd(s) 

t&j/ej-Fiji  l-ej)lejH  =  £(S& -F^aTf 
j=  i  j=  i 

Equation  of  the  transport  of  turbulent  kinetic  energy 
p V  •  uk  =  V  •  (p  +  pJak)Vk  +  Pk  +  Gk  —  pe 
Equation  of  the  turbulent  energy  dissipation 


(5) 

(6) 

(7) 


(cr 


pV  *  u£  =  V  •  (p  +  pt/(Te) Ve  +  ^C!  •  Pk  +  C\(\  —  C3)Gk  —  C2  *  pe 


1  + 


e/A: 


(8) 

in  which 

pt  =  pC^k2/s,  crt  =  ptCp/At,  Pk  =  pt(2Vw  •  Vw  +  (dw/ dy +  da/dx)2) 

Gk  =  gptpVT/(jt,  rj  =  ( Pk/gt)°-5k/s 

where  the  pressure  defect  p  is  the  difference  between  the  channel  and  the  ambient 
pressure  at  the  same  elevation,  and  Ta,  is  the  ambient  temperature  taken  as  the 
reference  temperature  for  this  study.  The  following  standard  constants  are  used  in 
the  k-s  model  using  the  Renormalization  Group,  RNG  method,  see  Yakot  et  al.[3] 

q,  =  0.0865,  Cj  =  1.45,  C2  =  1.83,  C3  =  0.8,  <rt  =  0.9, 

crk  =  0.8,  oE  =  1.15,  k  =  0.41,  ^  =  0.1695,  r]o  =  4.618 

Laminar  flows  can  be  derived  by  setting  pt  =  0  and  omitting  the  equations  for  k 
and  s.  The  associated  boundary  conditions  for  velocity  and  temperature  are  (see  Fig. 
1): 

at  the  inflow 

T=T.d  (9) 

on  the  insulated  wall 

u  =  v  =  0,  dT/dx  =  0  (10) 

on  the  heated  wall 

u  =  v  =  0,  q"  =  q"n  (11) 

At  the  first  inner  grid  point,  the  following  standard  equations  are  used  to  predict  the 
turbulent  kinetic  energy,  km  =  a  •  Ufn,  and  its  dissipation,  sin  =  pC^(k2/Rg).  The  fol¬ 
lowing  standard  constants  are  used  in  the  above  equations,  a  =  0.01  and  R  =  10. 

The  governing  equations,  eqns  1-8,  for  the  present  study  are  solved  by  the  finite 
element  code  FIDAP  with  a  pressure -projection  method  introduced  in  Fluid  Dynam¬ 
ics  International. [4]  The  segregated  iterative  solver  is  used  to  solve  the  transport 


B.  Moshfegh,  M.  Sandberg / Renewable  and  Sustainable  Energy  Reviews  2  ( 1998 )  287-301 


293 


variables  u,  v,  T,  k,  s,  and  p.  To  overcome  the  difficulty  of  numerical  convergence  and 
instability,  solutions  are  obtained  by  employing  hybrid  upwinding  and  relaxation 
schemes.  The  hybrid  relaxation  factors  are  maximum  0.1  and  minimum  0.01  for 
u,  v,  k,  and  s,  and  maximum  0.05  and  minimum  0.01  for,  T.  The  local  criterion  for 
numerical  convergence,  i.e.  the  maximum  relative  difference  between  two  consecutive 
iterations  for  any  local  variable,  is  less  than  10-3.  Convergence  is  typically  achieved 
after  about  250  iterations  with  a  central  processing  time  of  four  h  on  the  Dec  3000- 
900  AXP  computer. 

Non-uniform  grids  are  used  both  in  v  and  y-directions.  A  large  number  of  grid 
points  are  placed  at  geometrically  decreasing  distances  in  the  region  next  to  the  walls, 
where  specially  large  gradients  of  velocity  and  thermal  gradients  are  expected.  Four- 
node  quadrilateral  elements  are  used.  A  different  number  of  grid  points  are  tested  and 
the  resulting  flow  rate  and  heat  balance  from  these  tests  are  compared,  to  check  the 
accuracy  of  the  solution.  A  grid  system  of  79  x  231  nodes  is  found  sufficient  for  the 
solutions  presented  here. 


4.  Numerical  and  experimental  results 

The  numerical  results  presented  here  are  obtained  for  constant  uniform  input  heat 
flux,  q[ varying  from  20  to  500  W/m2.  Effects  of  the  emissivity  of  the  bounding 
surfaces  are  also  analysed  on  the  thermal  behaviour  of  the  air  passing  through  the 
gap.  The  ambient  air  temperature  is  293  K  and  the  Prandtl  number  is  Pr  =  0.708. 
Other  important  parameters  are  the  height  of  the  heated  wall  and  the  insulated  wall 
which  are  6.5  and  7.0  m,  respectively.  The  thermal  properties  of  the  air  are  calculated 
at  the  reference  temperature  of  300  K. 

4.1.  Flow  and  heat  transfer  at  different  constant  uniform  input  heat  flux 

The  outlet  air  velocity  and  outlet  air  temperature  profiles  as  well  as  the  surface 
temperatures  of  the  heated  and  insulated  walls  of  the  air  gap  subjected  to  a  uniform 
input  heat  flux  20,  300,  and  500  W/m2  are  presented  in  Figs  2  and  3.  The  corre¬ 
sponding,  Rad  =  g Pq"nd4 /vaHk,  for  the  present  cases  are,  2  x  1 0s,  3  x  109  and  5  x  109, 
respectively.  These  Rad  numbers  are  based  on  the  separation  distance  of  the  gap, 
d  =  0.23  m. 

Figure  2(a)  and  (b)  show  the  outlet  air  velocity  and  temperature  profiles.  At  20 
W/m2,  the  flow  is  assumed  to  be  laminar.  It  can  be  noted  that  the  convective  flow 
induced  by  the  horizontal  temperature  gradient  is  very  weak,  and  it  will  be  shown 
later  that  the  radiation  heat  transfer  between  surfaces  is  a  dominating  mechanism. 
Two  boundary  layers  are  developed  along  the  heated  and  insulated  walls.  At  a  large 
part  of  the  outlet  cross  section  area  the  outlet  air  temperature  is  equal  to  the  ambient 
temperature.  By  increasing  the  Rad  number,  the  flow  characteristic  of  the  air  will  be 
changed.  For  300  and  500  W/m2  the  flow  is  assumed  to  be  turbulent.  This  case  is 
dominated  mainly  by  the  turbulent  convection  heat  transfer.  This  effect  is  due  to  the 
fact  that  the  heat  is  mostly  convected  from  the  heated  wall  to  the  insulated  wall  by  a 
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Fig.  2.  (a)  p-component  velocity  profile  of  the  air  at  the  outlet  of  the  gap,  v0,  and  (b)  the  outlet  air 
temperature  profiles,  T0  —  Tr 
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Fig.  3.  Surface  temperature  of  the  heated  and  insulated  walls,  (a)  q"n  =  20  W/m2  and  (b)  q"n  =  300  and  500 
W/m2. 
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very  strong  turbulent  convection  along  the  horizontal  cross  section  area.  According 
to  Fig.  2(a),  there  is  two-peak  velocity  close  to  the  heated  and  insulated  walls,  but  the 
velocity  profile  of  turbulent  case  has  a  tendency  to  be  more  flattened. 

Figure  3(a)  shows  the  surface  temperature  of  the  heated  and  insulated  walls  for 
q"n  =  20  W/m2.  The  experimental  results  are  also  presented  in  Fig.  3(a).  The  instability 
tendencies  in  temperature  values  at  the  lower  part  of  the  heated  wall  can  be  a  reason 
of  the  separation  effect  close  to  the  lower  edge  of  the  heated  wall.  For  the  turbulent 
cases  (q"  =  300  and  500  W/m2)  Fig.  3(b),  the  surface  temperatures  increase  very  fast 
along  the  height  of  the  channel,  compared  to  the  laminar  case  where  the  changes  are 
very  slow.  It  is  also  interesting  to  note  that  the  temperature  difference,  between  the 
heated  and  insulated  wall  along  the  height  of  the  channel  is  almost  constant  for  the 
laminar  case,  while  for  the  turbulent  flow  the  temperature  differences  increase  along 
the  height  of  the  channel.  The  numerical  predictions  agree  well  with  the  results  from 
the  measurements  for  the  laminar  case,  as  shown  in  Fig.  3(a). 

Figure  4(a)  and  (b)  show  the  variation  of  the  average  outlet  air  velocity  and  the 
outlet  air  temperature  as  a  function  of  the  input  heat  flux.  The  outlet  air  velocity  and 
temperature  are  found  to  have  a  linear  relation  with  the  input  heat  flux  for  the 
turbulent  cases  in  the  log-log  scale.  A  linear  regression  analysis  is  performed  in  the 
log-log  scale  and  the  fitted  results  are  presented  in  the  legends  of  Fig.  4(a)  and  (b).  It 
can  be  seen  that  the  power  law  expression  provides  an  excellent  representation  of  the 
numerical  results.  For  comparison  the  experimental  results  and  the  lumped  parameter 
analysis  are  also  included  in  Fig.  4(a)  and  (b). 

Based  on  the  lumped  parameter  analysis  for  fully  turbulent  flow,  the  power  law 
relation  between  velocity  and  input  heat,  and  between  the  temperature  rise  at  the 
outlet  of  the  gap  and  the  input  heat  raised  to  an  exponent  equal  to  1/3  and  2/3 
respectively,  see  eqns  (12)  and  (13),  compared  to  0.353  and  0.654  for  the  numerical 
predictions.  A  more  detailed  analysis  of  the  lumped  parameter  method  will  soon  be 
presented  by  authors.  For  the  experimental  data  for,  q"n,  varying  from  200  to  300 
W/m2  the  exponent  is  0.34  for  the  relation  between  outlet  velocity  and  input  heat. 
This  is  within  the  measurement  accuracy  equal  to  the  theoretical  value. 


(12) 


A  72  =  Tn 


r.=-Lph_ 

90  la HA? 


1/3 


#2/3 


(13) 


For  convenience  the  following  terms  have  been  introduced,  specific  buoyancy 
flux,  B  =  g/3q/pCp,  and,  \Jj  =  p(l/p0  —  l/pd  +  XH/d+0.5(l  +  kentr).  The  result  from  the 
lumped  parameter  analysis  is  based  on  the  friction  coefficient,  X  =  0.024,  and  the 
entrance  loss,  kentr  =  0.5,  which  is  the  standard  value  for  a  sharp-edged  inlet.  These 
input  parameters  predict  a  velocity  somewhat  higher  than  the  recorded  velocities. 
One  plausible  explanation  is  that  the  entrance  losses  are  larger  because  the  inlet  is  a 
90°  bend,  see  Fig.  1. 

The  numerically  determined  average  heat  transfer  coefficient  on  the  heated  and 
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Fig.  4.  Comparison  of  experimentally  and  numerically  determined  outlet  air  velocity  (Fig.  4(a))  and 
temperature  (Fig.  4(b))  results.  The  square  shown  in  Fig.  4(b)  is  the  recorded  temperature  dilference  while 
the  filled  circle  shows  the  obtained  temperature  dilference  based  on  input  heat  and  recorded  flow  rate. 
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insulated  walls  at  a  different  input  heat  flux  are  plotted  in  Fig.  5  as  a  function  of,  q'f 
For  fully  turbulent  flow  cases,  i.e.  for,  q"n,  varying  from  200  to  500  W/m2  (cor¬ 
responding  to  a  RaH  number  varying  between  1.3  x  1015  to  3.2  x  1015  based  on  the 
height  of  the  heated  wall),  the  best  fitted  lines  showed  that  the  power  exponent 
between,  hh  and  q'f  and,  hc  and  q"n,  is  equal  to  0.247  and  0.228,  respectively.  The 
exponents  are  in  good  agreement  with  the  result  presented  in  Vliet  and  Liu. [5]  The 
results  are  presented  in  the  legends  of  Fig.  5. 

The  percentage  of  the  input  heat  which  is  dissipated  to  the  air  via  convection  and 
via  radiation  to  the  bounding  surfaces  is  show  in  Fig.  6.  The  results  show  that  the 
magnitude  of  the  radiative  heat  transfer  is  comparable  to  the  convective  term  at  low, 
q"n,  but  it  has  a  tendency  to  decrease  by  increasing  the  input  heat  flux.  For  input  heat 
flux  200  W/m2  more  than  67%  of  the  input  heat  is  transferred  to  the  air  via  convection 
and  it  is  increasing  very  slowly  by  increased  input  heat  flux. 

4.2.  Flow  and  heat  transfer  at  different  emissivity  of  the  bounding  surfaces 

The  relative  flow  rate  through  the  gap  decreases  by  reducing  the  emissivity  of  the 
bounding  surfaces,  but  on  the  other  hand,  the  average  outlet  air  temperature  increases 
(the  input  heat  flux  is  constant),  see  Table  1.  In  other  words,  the  cooling  capacity  of 
the  cavity  will  be  reduced  and  result  in  a  lower  efficiency  for  the  panel.  It  can  be  seen 
also  from  Table  1  that  the  emissivity  of  the  bounding  surfaces  strongly  influences  the 
surface  temperature  of  the  insulated  and  heated  walls.  The  difference  between  the 
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Fig.  5.  Average  heat  transfer  coefficient  for  the  heated  and  insulated  wall. 
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Fig.  6.  The  amount  of  the  heat  transferred  from  the  heated  wall  to  the  air  by  means  of  convection  and 
radiation. 
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average  surface  temperature  of  the  heated  wall  for  case  one  and  four  is  more  than 
14°C.  This  temperature  rise  will  reduce  the  efficiency  of  the  solar  cell  by  0.6%.  The 
fact  is  that  the  amount  of  radiative  heat  transferred  from  the  heated  to  the  insulated 
wall  will  be  reduced  by  increasing  the  thermal  resistance  due  to  the  radiation.  For 
surfaces  with  a  lower  emissivity  factor  most  of  the  input  heat  will  be  transferred 
directly  from  the  heated  wall  via  convection  to  the  air.  For  the  pure  convection  case 
(eh  =  sc  =  0)  the  radiative  heat  flux  is  zero  and  the  consequence  is  that  the  convection 
term  for  the  insulated  wall  can  be  neglected,  the  surface  temperature  of  the  insulated 
wall  is  close  to  the  air  temperature  adjacent  to  the  wall  and  the  surface  temperature 
of  the  heated  wall  will  rise  to  the  highest  possible  temperature. 

The  results  also  indicate  that  the  amount  of  radiation  heat  penetration  through  gap 
openings  is  also  reduced  from  case  one  to  four.  In  other  words,  using  surfaces  with 
lower  emissivity  factor  will  allow  the  convective  heat  transfer  to  dominate  the  heat 
transfer  mechanism  while  the  radiation  heat  exchange  between  surface  boundaries 
will  be  diminished.  The  effect  of  the  separation  distance,  d,  on  the  thermal  behaviour 
of  the  air  was  presented  in  an  earlier  report  by  the  authors,  see  Moshfegh  and 
Sandberg[6]  and  Sandberg  and  Moshfegh[7j. 

5.  Conclusions 

An  investigation  was  made  to  provide  better  understanding  of  the  mechanism  of 
heat  transfer  in  the  air  gap  behind  PV  panels  by  means  of  both  numerical  and 


Table  1 

Effect  of  emissivity  of  the  bounding  surfaces  on  the  channel  flow  and  thermal  characteristics 
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{do  ddconvl  din 
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(do+ddrJdm 
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1.0 

1.0 

27.8 

0.89 

50.6 

35.0 

6.7 

5.8 

68.0 

32.0 

97.8 

2.2 

2 

1.0 

0.1 

28.3 

0.83 

62.7 

24.3 

6.3 

5.0 

89.4 

10.6 

96.6 

3.4 

3 

0.1 

1.0 

28.5 

0.84 

64.2 

23.9 

6.3 

5.6 

93.4 

6.6 

99.8 

0.2 

4 

0.1 

0.1 

28.7 

0.82 

65.2 

22.0 

6.4 

4.9 

96.2 

3.8 

99.8 

0.2 

d  =  0.23  m ,H=  6.5  m,  Pr  =  0.708,  Rad  =  3.109,  T{  =  Ta  =  20°C,  q"  =  300  W/m2 
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experimental  method.  For  vertical  channels,  UQ  and  Ta  increases  by  increasing  the 
input  heat  flux.  However,  for  a  fixed  input  heat  flux,  UQ  decreases  by  using  surfaces 
with  lower  emissivity  factor,  due  to  the  increments  in  temperatures  at  the  outlet  of 
the  channel,  T0.  In  terms  of  efficiency  of  the  solar  panel,  an  indication  that  is  beneficial 
to  have  surfaces  with  high  emissivity  factor. 

For  fully  turbulent  flow,  UQ  and  T0,  are  according  to  the  numerical  predictions 
proportional  to  the,  q"nm,  where  m  is  0.353  and  0.654,  respectively.  For  the 
200  ^  q"n  ^  500  W/m2,  hh  and  hcccq''Qn,  where  n  is  0.247  and  0.228,  respectively. 
According  to  the  measurements  for  the  interval  200  <  q"n  ^  300  W/m2  the  exponent, 
m,  is  0.34  for  UQ.  The  exponents  show  good  agreement  between  the  numerical  pre¬ 
dicted  results  and  the  experimental  results. 

The  analysis  revealed  the  importance  of  radiation  heat  exchanges  for  the  heat 
transfer  mechanisms  between  the  channel  walls.  For  input  heat  flux  equal  or  more 
than  200  W/m2  nearly  30%  of  the  input  heat  is  transferred  to  the  otherwise  unheated 
wall  via  radiation  and  then  it  is  dissipated  to  the  air  through  convection. 
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